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SUMMARY 


The  reliability  theory  for  fatigue  lives  of  one-crack,  two-stage  models  with  hijacking 
has  been  extended  to  allow  for  arbitrary  inspection  intervals  with  possible  arbitrary 
renewals.  The  life  distribution  is  similar  with  a  modified  initiation  density  and  the  associated 
moment  generating  function  is  formally  identical.  It  is  possible  to  identify  an  unsteady 
Markov  chain  formed  by  the  combined  fatigue  and  inspection  process. 

Two  FORTRAN  programs  have  been  developed  from  this  formulation,  one  of  which 
allows  for  random  crack  rates  and  run-time  setting  of  inspections.  w 
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ABSTRACT 

The  reliability  theory  for  fatigue  lives  of  one-crack ,  two-stage  models  with  hijacking 
has  been  extended  to  allow  for  arbitrary  inspection  intervals  with  possible  arbitrary 
renewals.  The  life  distribution  is  similar  with  a  modified  initiation  density  and  the  associ¬ 
ated  moment  generating  function  is  formally  identical.  It  is  possible  to  identify  an  unsteady 
Markov  chain  formed  by  the  combined  fatigue  and  inspection  process. 

Two  FORTRAN  programs  have  been  developed  from  this  formulation,  one  of  which 
allows  for  random  crack  rates  and  run-time  setting  of  inspections. 
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DISTRIBUTION 


NOTATION 


a  (a) 

do 

B 

Bi  =  [M 
£ 

f(a\t) 

fa{U) 

m 

MO 

£ 

Fa 

Fi 

gi,  gt(t  -  Tt) 

Gt 

h^t0)  =  e'oiof^t,,) 

Hi 

I,J,K 

M 

Mi 

Nt 

P,  Pi(a)  =  />,('.) 

Pi  =  P(Tt  -  T,o) 

Pio 

Pi  =  PAT,  -/  +  /.) 
P)=  1  -P/Pra 
Pr(.  |.) 

PAu\to) 


Crack  size(s). 

Initial  crack  size. 

Infinite  transition  matrix,  equation  (4.12B). 

Fifth  order  transition  matrix  for  the  inspection  at  time  7|. 
Expectation. 

Conditional  probability  density  of  crack  length 
(see  <£(.),  0(.))  at  time  t. 

Reliability  type  density  of  crack  life  (see  Equation  (2.6)). 

Density  of  initial  life  in  the  absence  of  any  inspection  or  repair. 
Initial  life  density  as  affected  by  inspection  or  repair. 
fj(t)  defined  in  (Tt,  7i+i]  extended  to  (Tt,  oo). 

Cumulative  density  of  distribution  function,  not  necessarily  proper. 

1  -  Fa. 

1  - 

Density  of  initial  lives  in  structures  repaired  after  inspection  Tt. 

(This  is  initial  data). 

Distribution  corresponding  to  gt,  not  necessarily  proper. 

Used  for  Equation  (5. ISA). 

Defined  by  (S.10A).  This  is  a  hypothetical  distribution  of  initial  lives 
from  a  hijacked  distribution  for  structures  repaired  at  Tt. 

Terms  of  (4.7);  see  (4.10). 

Generic  form  of  moment  generating  function,  E.  exp(— u). 

Moment  generating  functions  from  g See  Section  4.4. 

Defined  by  (5.15)  and  (5.15A). 

Probability  of  rejection  for  an  inspected  structure  with  a  crack  of 
length  a  or  at  time  ta  after  initiation  (operating  characteristic). 
Normally  used  without  subscripts. 

Used  in  Section  4.3. 

Laplace  transform  of  P.  with  respect  to  ta  when  t0  =  Tu,  +  i.e.  when 
the  crack  starts  just  after  an  inspection. 

The  probability  of  rejection  at  an  intermediate  inspection  (Section  3.2) 
Probability,  absolute  or  conditional. 

Defined  by  (4.2).  Overall  acceptance  probability  during  growth  of 
some  crack. 


Qi  Abbreviation  used  in  /,  Equation  (4.8)  and  preceding. 

R  =  R(a)  =  R(im)  Average  local  crack  growth  rate. 

Ri  Probability  of  rejection  at  i  th  inspection,  time  7< .  Equation  (3.1). 

R  =  {/ti|<  =  1, ...,  oo}  Infinite  vector  of  &'$. 

r(a)  =  r(/«)  Risk  function. 

r0  =  r(0)  Hijack  risk  component. 

£(/)  Arbitrary  step  function,  Section  4.2. 

S(u)  Laplace  transform  of  5(f). 

t  Current  equivalent  time  or  cycles. 

/.  Other  times  according  to  subscript 

7}  Known  inspection  times,  T„  =  0. 

u,  v,  m>  Laplace  transform  variables. 

ATt  =  7|+i  —  Tt  i-th  inspection  interval. 

y,  8  Real  shifts  for  Bromwich  inversion  contours. 

<Kt)  Distribution  of  final  lives  without  inspection,  (2.6). 

MO  Distribution  of  Anal  lives  with  inspection  and  repair,  (3.6). 

•  Convolution  function. 


Subscripts 

The  subscripts  associate  the  main  symbol  with  the  quantity  listed. 
o  Crack  initiation  without  renewals. 

*  Crack  initiation  affected  by  renewals. 

•<  As  above  but  distinguished  by  inspection  period. 

a  Crack  growth  time. 

c  Set  complementation. 

i  Inspection  period  (7),  7}+iJ. 

f'o,  io  The  inspection  just  before  the  start  of  the  current  crack  i.e. 

Tte  <  to  ^  f|»+i 


Sets 

Union  and  intersection  are  here  denoted  by  “+”  and  the  normal  product 
convention  while  complements  are  indicated  by  the  subscript  c.  Figures  I  and  S 
illustrate  some  of  the  events. 

A  Attrition. 

E  Universe. 

M  =  CAcRc  Mainstream  —  cracked  but  still  in  use. 

R  Rejected  at  inspection, 

d  Empty  set  (see  Section  5.1). 


1.  INTRODUCTION 

It  is  a  truism  that  fatigue  life,  especially  with  single  cracks,  consists  of  initiation  time  and 
crick  growth  time.  The  models  based  on  this  approach  when  the  initiation  is  random,  have 
been  described  in  previous  Reports12  3.  In  the  last  of  these,  considerations  of  continuity  of 
probability,  together  with  deterministic  cracking,  led  to  a  first  order  partial  differential  equation 
(true  for  vector  cracks  also),  which  is  the  same  as  the  continuity  condition  for  compressible 
flow.  It  is  also  a  degenerate  form  of  the  second  order  Fokker-Planck  equation  which  it  would 
become  for  random  cracking.  Most  generally  this  describes  the  infinitesimal  evolution  of 
probabilities  associated  with  continuous  Markov  processes4'5. 

In  an  earlier  Report3  the  density  of  crack  lengths  (and  thence  that  of  total  life)  were  found 
from  the  continuity  equation  without  considering  inspections,  though  these  were  mentioned 
briefly.  Under  the  term  hijacking  it  also  introduced  the  effect  of  losses  not  due  to  fatigue.  The 
present  Report  extends  the  previous  solution  to  include  inspections  as  well  as  the  hijack  risk. 
It  considers  the  distribution  of  total  life,  the  moment  generating  function,  and  the  transition 
matrix  for  changes  of  state  between  inspections. 


2.  PREVIOUS  RESULTS 

Before  proceeding,  we  shall  summarise  results  for  the  one-crack  model  without  inspection. 
As  before,  we  shall  use  the  generic  notation  /,  F,  <f>,  and  M ,  to  denote  density,  distribution 
functions  and  moment  generating  functions  of  their  arguments.  Density  is  affected  by  attrition 
and  where  necessary,  these  symbols  will  be  subscripted  to  avoid  confusion. 

For  vector  cracks  the  continuity  equation 

^=/(div  R(a)  +  r(a))  (2.1) 


holds,  where  D  denotes  total  derivative  and  R(a)  =  d»/dt,  the  previously  averaged  crack  rate, 

r(a)  =  total  risk  function  including  hijack  risk,  and 

div  R(a)  =  dRjda  for  a  single  cracking,  a  known  function  of  crack  length. 

When  this  is  expanded,  one  obtains  the  degenerate  Fokker-Planck  equation  (for  single 
cracks) 

<12) 

with  the  characteristic  equations 

dt  =  da/R(a)  =  -df/f(dR/da  +  r(a)).  (2.3) 

In  this  the  crack  trajectories  are  characteristics  and  the  general  solution  for  crack  length 
density  is 


f — f  (°o\to) 


R(ao) 

R(a) 


—  f(do\to) 


*(«.) 

R(o) 


(2.4) 


introducing  the  growth  time  /■  to  reach  crack  length  a.  Here  at  is  the  initial  crack  size,  a  constant, 
and  t0,  the  initial  life,  corresponds  to  the  crack  a\t.  The  growth  time 
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t  o  —  t  —  to  • 

When  the  hijack  risks  are  denoted  r0>  and  the  boundary  conditions  are  included  r(t)  =  r(a(t)) 

e~r°tofo{to 


f(a\t)  = 


m 


^exp|-|\(/Vtj 


where  the  generic  /  refers  to  different  density  functions  according  to  its  subscript. 
It  is  now  possible  to  average  the  risk  r(a )  at  time  t  for  the  overall  life  density. 

<K0  =  roe  r°to  [1  -  F0(/)]  +  e~Tot°fo(to)  */«(/«) 

'ta 

where /«(/«)  =  /■(/„)  exp  —  I  r(t)dt,  the  reliability  based  crack  life  density. 


j: 


The  moment  generating  function,  exp  —  ut,  follows  as 


M^u)  =  Mo(u  +  r0)  Ma(u)  + 


— r 


1  —  Mo(u  +  r0) 


(2.5) 


(2.6) 


(2.7) 


which  differs  slightly  from  the  corresponding  equation  of  Reference  3  since  —u  is  used  here 
to  assure  convergence  for  positive  u. 


2.1  Effect  of  Inspections 

This  is  twofold:  in  the  first  place,  the  density  of  crack  length  will  be  reduced  suddenly 
at  a  number  of  steps  corresponding  to  each  inspection.  Secondly,  structures  rejected  at  an 
inspection  may  be  returned,  after  repair,  to  the  population.  This  brings  the  whole  problem 
into  a  close  relationship  with  statistical  renewal  theory.  However,  the  latter  does  not  include 
hijack  risks,  nor  the  two  stages  involved  in  fatigue. 

We  will  first  consider  renewals  and  then  the  corresponding  moment  generating  functions. 
In  the  following,  it  is  most  convenient  to  have  all  variables  as  time,  and  continue  to  use  t0,  ta 
and  t  to  denote  initiation  time,  crack  period  and  their  sum.  The  subscript  *  will  refer  to  densities 
or  generating  functions  affected  by  renewals.  Inspection  times  are  7),  i  =  0,  1,2  etc.  where 
T0  =  0.  The  i-th  inspection  interval  and  the  associated  quantities  occur  after  7). 


Repairs 


Not  cracked 


FIG  1.  POSSIBLE  STATES  OF  STRUCTURES 
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3.  LIFE  DISTRIBUTION 

At  a  particular  time  the  population  of  structures  suffering  two-stage  fatigue  with  inspections 
and  attrition  may  be  described  by  the  Venn  diagram  of  Figure  I.  In  the  course  of  time  the 
arrows  indicate  the  evolution  of  members  of  the  population;  one  may  equivalently  imagine 
uniform  measure  density  over  the  rectangle  and  the  subset  boundaries  moving  left  and  inwards 
opposite  to  the  arrows.  The  zero  set  corresponds  to  rejection  of  uncracked  structures.  This  is 
the  general  process  we  intend  to  describe. 


3.1  Renewal  Allowance  for  Single  Cracks 

At  Ti  —  just  before  inspection,  suppose  the  partial  density  of  vector  crack  lengths  is 
/(a|7(— ),  which  therefore  becomes 

/(a|7i  +  )  =  (l  —  /><(a))/(aiTi— ) 

just  afterwards.  With  respect  to  the  whole  population  this  gives  the  probability  of  rejection  at 
the  i-th  inspection  as 

f 

Pi(*)f(*\Tt-)d*  (3.1) 

JO 

We  also  define  Pi  the  operating  characteristic  of  the  inspection  method,  in  terms  of  growth  time 
so  that  Pt(a )  =  />((fo).  The  context  will  indicate  which  definition  is  current. 

After  Tt,  the  population  will  also  include  repaired  or  modified  structures  returned  from 
inspection.  It  may  also  be  depleted  by  those  retired  from  service.  Statistically,  these  repairs  or 
retirement  affect  the  initial  life  distribution.  After  each  inspection  we  will  assume  that  structures 
are  repaired  to  the  same  condition  (not  necessarily  the  original  as  new  condition)  despite  the 
possibility  of  differing  crack  lengths  being  discovered.  (Since  it  is  mathematically  easy,  the 
restored  condition  shall  be  assumed  specific  to  the  particular  inspection). 

This  assumption  requires  us  to  first  consider  the  density  of  crack  initiation  which  may 
determine  boundary  conditions  for  f(a\t  >  Ti)  (see  Ref.  3). 

For  one  crack  in  the  absence  of  inspections,  let  e  Tot(\  —  Fa(l))  be  the  survivorship  function 
for  initial  life  when  hijack  attrition  is  included  (Ref.  3).  When  inspections  are  included  but  not, 
for  the  moment,  hijacking,  let  the  corresponding  survivorship  fraction  (based  on  the  original 
population)  be  F.<_i(r)  for  t  <  Tt  say. 

Just  after  an  inspection  all  the  uncracked  structures  will  be  retained  and  also  augmented 
by  the  fraction  J?t(l  —  C<(0)]  of  inspected  and  repaired  structures. 

Symbolically 

1  -  F.t  (7)+)  =  [I  -  F.t_i  (Tt-)]  +  R,[l  -  G,(0)]  (3.2) 

where  Gt(l)  is  a  general,  possibly  partial,  distribution  function  of  initiation  at  times  T<+  t  of 
structures  repaired  at  time  Tt.  At  time  t,  Tt  <  t  <  T(+ 1  with  hijacking  try 

e-rot  (l  _  F.,(t))  =  e-r<>i  (j  _  F.,_i(/))  +  Rte-'o»-Tn  (|  _  ct(t  -  T()) 


or 

Ft  =  Ft. ,  ( t )  +  Hi 

where  F  —  1  —  F  and 

Hi  =  er°Tl  R,[l  -  Gi(t  -  7,)] 
Recursive  substitution  then  leads  to 


/ 

Ft(t)  =  F0  +  '^H, 

j- 1 


which  requires  an  interpretation  of  F.0,  G0  etc. 


(3.3) 


(3.3A) 
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Initially  (3.3)  indicates  that 

1  -  F.„(T»)  =  1  -  F.  xiTo)  +  e’-o'1'-  *„[i  -  G„(0)]  (3.3B) 

and  we  would  expect  that 

A0(/|/  <  Ti)  =  F„U), 

the  initial  failure  distribution.  If  we  define  T0  =  0  then  for  0  <  /  <  7i,  (3.3)  becomes 

I  -  F.o{t)  =  I  -  F0(t)  =-  I  -  F.  i(t)  +  R0{\  -  G0(t  -  To)]  (3.3 C) 

The  definitions  may  now  be  arbitrary.  The  most  convenient  for  us  is  F._i(l)  =  Fd(t)  and 
R0  =  0  leaving  G0  arbitrary. 

(Note  that  all  F.u  Gt  etc.  are  defined  on  (0,  oo);  F.i  is  applied  on  (Ti,  7t+i)  and  F.  =  F.t  etc. 
everywhere  with  the  index  suited  to  t.)  To  summarise  T0,  R0  =  0,  F._i(f)  =  F.^t)  =  Foil) 
and  Go  is  arbitrary. 

In  (3.3)  or  (3.3A)  the  restored  factor  exp— r0t  allows  for  hijacking.  These  equations  now 
provide  some  of  the  boundary  conditions  for  f(a\Tt  <  t  <  Ti+i)  when  the  solution  has  marched 
to  Tt+i  —  .  Then  (3.1)  provides  allowing  (3.3C)  to  operate  over  the  next  inspection  inter¬ 
val  (7m»  Ttri). 


3.2  Density  Function  for  Attrition 

Now  include  the  hijack  factor  in  (3.3C)  and  consider  the  density  of  initiation 
d 

—  e-  r0t  f(  =  r0e  rul  Ft  +  e~rotf.i(t)  (3.4A) 

dt 


where  the  first  term  is  obviously  the  local  incidence  of  hijacking  and  therefore  a  direct  component 
of  the  density  <f>(t )  of  total  life. 

From  (3.3A)  the  second  term  is 


e  rot 


i 

2 

i-  l 


er°Ti  R)g)(t  -  T}) 


(3.4B) 


where  gj  is  the  initiation  density  corresponding  to  Gj  which  may  be  partial  or  defective. 

Consider  f(a\Tt  <  t  <  7j+i)  after  the  i-th  inspection.  Such  a  crack  would  have  begun  at 
t  —  t„( a)  and  would  therefore  have  been  inspected  an  integral  number  of  times.  The  operating 
characteristic  P(a)  of  any  crack  is  the  defined  in  terms  of  crack  length.  In  keeping  with  our 
present  approach,  it  is  more  convenient  to  define  this  in  terms  of  the  growth  time  ta  as  P(ta). 
We  further  define  P]  as  the  characteristic  corresponding  to  the  y'-th  inspection  of  a  crack. 


FIG  2.  NOTATION  FOR  INSPECTION  DURING  CRACK  GROWTH 
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In  the  notation  of  the  figure 


Pi  =  P(T}  -  t  +  /.) . 


For  a  long  crack,  several  non-trivial  Pi's  are  possible.  Unless  specifically  stated,  several  Pj's  in 
the  same  expression  refer  to  the  same  crack  trajectory. 

Now  consider  the  basic  density  f(a\t)  allowing  attrition  of  cracked  structures  and  also 
inspections.  If  7<  <  t  <  7<+i  integration  (Ref.  3)  along  the  characteristic  and  multiplications  by 
Pj  =  1  —  P  produce 


R(a)f(*\t)  =  Pt0+ 1 ...  Pt  exp 


e  r°t°f^0(t  -  U) 


(D'M-M'il 


roTJRlg)(to-  T,)  (3.5) 


where  ta  is  the  growth  time  and  Ti0<  t  —  ta<  Ti0+i .  As  before,  with  the  hijacking,  this  produces 
the  attrition  density 


t  <  Ti+i)  =  roe  r°‘Fi  +  |  r(t)f(%\t)da 

do 


i. 


oo 


=  roe-r<>t  [ i  -F„(r)]  + 


f  e-'o'<Ut0)  FT  F(Tj 

J  0  ,.+> 


to)dFa(ta) 


(3.6) 


with  dFtt(ta)  from  (2.6),  the  initiation  t0=  t  —  ta,  as  always,  and  Pj  is  described  in  full.  In  (3.5) 
the  finite  upper  limit  t  is  no  restriction  because  /.=  0  for  negative  arguments.  Recognition  of 
this  will  aid  the  manipulations  below. 


4,  MOMENT  GENERATING  FUNCTIONS 


Equation  (3.6)  appears  to  be  a  convolution  but  the  status  of  Pj  is  uncertain. 
Let  us  form  the  moment  generating  function  as 


z 

1  =  0 


'T,+i 

<f>  (r|i)  e  uldt 

J  Tt 


r0e-'ro+«)t  [i  _  F.(t)]dr  f 


^  CTt+i 


y  r^.)«»-«^(9)n  P(Ti-l„)dFa(ta)dl, 

tlJTi  Jo  \+r 

where /.  is  the  global  form  of  /.(,  the  whole  of  the  modified  initial  life  distribution. 


(4.1) 


4.1  Inspection  Functions 

Consider  the  product  of  the  inspection  factors  as  functions  of  initiation  l0  and  of  growth 
time  la-  For  fixed  t0  this  term  is  a  step  function  of  ta  (determining  i  —  i0)  but  it  is  at  least  piecewise 
continuous  in  t0  since,  almost  always,  7)—  t0  continuously  determines  crack  sizes  at  inspection. 

Occasionally,  through  i0  decreasing,  t0  introduces  another  factor,  but  if  Pi(aB)  =  0,  which 
we  now  assume,  continuity  with  respect  to  t0  will  be  retained. 

Let  us  abbreviate  the  inspection  factor 

i 

]~J  (I  —  P(Tj—  to))  =  Pdia It.)  =  P.  (4.2) 

/.+  1 

where  P.  is  a  step  function  decreasing  from  unity.  We  now  intend  to  treat  (4.1)  as  a  convolution 
by  taking  some  of  the  inspection  terms  as  part  of  the  crack  life  distribution.  In  (4.1)  the  region 
of  integration  is  the  infinite  sector  shown  in  Figure  3.  The  operation  began  with  integration 
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along  strips  such  as  AB  which  were  then  extended  to  horizontal  inspection  bands  whose  contri 
butions  were  finally  summed. 


FIG.  3  REGION  OF  INTEGRATION  FOR  MOMENT  GENERATING  FUNCTION 

If  we  change  the  variable  t  to  l0 ,  then  we  may  do  the  first  integration  along  CD  and  include 
the  summation  by  adopting  the  limits  of  (0,  co)  for  ta- 
Thequus  with  (4.2)  ation  (4.1)  becomes 

M£u)  =  (1  -  M.(u  +  r0)) 

u+  r0 

Tf+i  f"  co 

f'(t0)e~ lu  ‘  r°)to  P,e  uta  dFa{ta)dt0  (4.3) 

Ti  Jo 

where  M.=  E,  exp  (—ut0)  =  MGF  of/.. 

When  P.  =  1  the  integral  reduces  to  the  uninspected  form  M.(u  +  r0)Ma(u)\  in  general, 
these  terms  introduce  convolutions  of  transforms. 

4.2  Step  Function  Transforms 

We  have  already  seen  that  P.  is  a  step  function  for  given  t0  and  in  the  figure  for  (4.1)  the 
first  integral  along  CD  traverses  an  infinite  number  of  inspections.  In  (4.3)  then  P.,  given  by 
(4.2)  is  interpreted  for  all  Thus  P.(ia\i0)fa(ta)  may  be  regarded  as  a  defective  conditional 
density  of  cracking  life. 

By  the  convolution  theorem  for  Laplace  transforms 

o°  |  fy  +  i°° 

P.(ta)fa(ta)e  ui°dta  =  ^.  P.(u  -  v)Ma(u)dv  (4.4) 

0  27,1  J  y  -  »oo 

where  the  bar  indicates  a  Laplace  transform,  as  is  Ma- 
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For  any  step  function  S(t)  such  as  that  below 


FIG  4.  TYPICAL  STEP  FUNCTION 


it  is  easily  shown  that,  in  the  notation  indicated, 

00 

£(«)  =  ul  2  (A  -  Pi-i)e  ut‘ 

i=  o 

In  the  present  instance  with  steps  from  (4.2),  one  would  have 

1 1  =  Tto -fl  to,  ...  ,  tj  =  Tio+J  to. 


(4.5) 


with  t0  arbitrary.  Part  of  the  t0  dependence  is  in  the  value  of  7)o+i  —  to- 

Let  P(0  be  the  transform  of  P.  corresponding  to  t0  =  Ti0  +  c.  Then  for  r0  elsewhere  in 
this  inspection  interval,  the  transform  of  the  step  function  is 

e“l‘o-Ti0>P(0  +  M-l(l  _  euito-Tio^  (4.6) 

The  first  term  has  a  simple  shift  operator  but  the  second  “spillage”  term  is  required  because 
the  backwards  shift  truncates  S(t )  for  t  <  0. 


4.3  Complete  Moment  Generating  Function 

We  now  perform  the  transforms  indicated  by  (4.3)  requiring  the  transform  of  step  functions. 
From  (4.5) 


oo 

Pio  =  M-1  (Pi—  P(-l)e~u{Tl~Ti0\  Pi  =  1,  Pto  =  1,  Pio- 1  =  0, 


giving  (4.6)  as 


with 


■t o  u  i  2  (Pi-  Pi- i)<rwr°  +  w-Hl  -  eul°~uTi°) 

io 

Pi  =  P(T(-  Ti0). 


The  inner  integral  (4.4)  of  (4.3)  is  now 

!  Cv  +  «°° 

2 Hi 


1 

+  2ni 


y-ievto  (P(  —  P(_i)e  vTi  Ma(u  —  v)dv 
Y  ~  '00  ~ 

‘'y  +  icc  j  1  fy+ioo 


Ma(v)dv  -  . 

U  —  v  2ni 


y  —  ico 

with  the  u  —  v  argument  variously  placed. 


J 


y-lgViTo-Tlo'  M  a(u  —  v)dv  (4.7) 


y  —  ICO 
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All  the  components  of  (4.7)  are  similar  so  that  we  now  substitute  the  first  into  (4.3)  to  obtain 


OD 

where  Qi0  =  ^  (Pi~  Pi-i)e~vTllv. 


f4to)e~{,l~v+roH°Ma(u  —  v)Qi0  dvdt0 


If  the  order  of  integration  is  reversed  then  in  each  interval  the  inspection  factor  is  constant 
with  respect  to  t0.  Thus  the  t„  limits  of  /  may  be  infinite  provided  the  inspection  factor  is  regarded 
as  a  step  function.  The  t0  integral  is  thus  the  transform  of  the  products  of /.  and  this  step  function. 
This  introduces  another  convolution  whence 


e+r0>T( 

(Qi  ~  Qi- 1) - : - dvdw 

u—v—w+r0 

1=0 


where  Q_i  —  0. 

In  the  w-plane  there  is  a  pole  at  u  —  v  +  r0.  Contour  integration  then  leaves  the  cancelling 
residues, 


7  =  2w'  Ma^U  ~  _  v  +  r»)  ^ 

V  1—0 


dv 

2<-0  7  "  0 


The  next  component  of  (4.7)  leads  to 


%Tu>+  i  j  *y  +  ice 
/  ir~.  f'(to)e 

t^0  J  T*o  J  Y-i° ° 


M+r»,,»A/0(M  —  v)  —  dt0 


=  s— .  A  fe(u  +  r0) 
2  m 


‘y  +  ioo 


y  —  ico 


Ma(u  —  v)  —  ,  absorbing  the  summation, 


=  Mdu  +  r0)HU)fa(t)  =  M.(u  +  r0)Ma(u)  (4.91 

using  the  convolution  theorem. 

In  the  third  component  of  (4.7)  the  presence  of  Tto  indicates  the  presence  of  another  step 
function.  In  (4.3)  it  leads  to 


K=  - 


Cy  +  i 

ri  I 

J  Y~  i' 


/.(fo) - Ma(u  —  v)e~vT‘°  dv 

—  ico  v 


J  Tio  J  Y  ~  ioo 

where  the  step  function  heights  are  v_1  exp  —vT<0.  In  the  same  way  as  before,  using  (4.5), 


1  f*  +  100  i  jV+io O  ^  e-tv-ic>T/o  —  e-<v-ui)Tio~l 

L  .  2ni  I  +  - TT7 - 

J  8  —  ico  J  y  —  ico  TTo 

Once  again  the  only  w-pole  is  for  w  =  v  in  which  case  the  residues  again  cancel;  K  =  0. 
When  all  these  results  are  assembled  ((4.3),  4.8)  and  (4.9)) 

Mj(u)  =  Hijack  term  +  /  +  J  +  K 


-7  fi- 

u  +  r0|_ 


M0(u  +  r0)  +  M.(u  -f  r0)Ma(u) 


which  is  exactly  the  same  as  it  would  be  without  inspections  except  that  M.  =  /.  and  includes 
renewals.  This  last  fact  means  in  general  that  3/^(0)  /  1 .  Other  cumulants  of  the  total  life  are 
shown  in  Table  1. 


TABLE  1 

Cumulants  of  Fatigue  Life  in  Terms  of  Initiation  and  Crack  Time  Moments 


M.(u)  =  MGF  of  cycles  to  initiation 


M.<*>(-r0)  =  +  C Kr0) 


where  /it  —  Arth  moment,  /i0  =  1  — 


00 


/tell  -  G*(ao)] 


Af«(w)  =  1  +  *1 U  ^a2U2  .  . 


*o  =  1, 


=  MGF  of  cracking  time,  moments  a*. 


Put 


where 


mt  =  M.ik)(r0)-,  At  =  at  —  kat-i/r0. 


r0  =  Risk  of  loads  above  ultimate  and  of  hijacking. 


Then  for  total  life  we  have  the  cumulants: 

xi  =  m0Ai  +  r0->;  (Mean) 

*2  =  2m\A\  +  m0Ai  —  m02A\2  4-  r0~2  (Variance) 

*3  -  3m2A\  +  3miA2  +  m0A3  —  6m,/niAi2  —  3m0zA\Ai  +  2m03Ai3  +  2/r03 
*4  =  4m3Ai  +■  6m2A2  +  4miAs  4-  m0^4  —  \2ntom2Ai2  —  24m</niA\A2  +  24m02miAi3 
-  \2nn2Ai2  -  m02{4AlA3  +  3 A22}  +  l2m03Ai2A2  -  6m0Mi4  +  6/r„* 


4.4  Transform  of  f. 

Though  informative,  equation  (4.10)  is  still  not  presented  in  terms  of  basic  MGFs,  relying 
as  it  has  on  M..  To  relate  M,  to  M0  and  other  known  MGFs  define  /.  piecewise  from  (3.4). 
Then  its  transform 


a  rrM 

M.(u)  =  2_j  e~u‘f.((l)dt 

i=o  J 

A  rrl+i  ' 

=  Mo(u)  4  £  >  eroTr«  R)g,(t  -  T,)dt 

(=o  J  Ti  | 

in  which  A/o(w)  -  E0  exp  —  ut0-  The  Rf,  the  overall  probabilities  of  rejection  at  times  7),  are 
effectively  constants. 

Changing  the  order  of  summation  absorbs  one  sum  to  produce 


\Uu)  =  MJu)  + 


*  r  °° 

V  RieroTJ  e~ut  gj(t  —  Tj)dl 

i  J  7) 


2- 

o 


e-fu-r„)Tj  A/;(W), 


defining  M/(u)  =  MGF  of  gj('). 


R0  =  I  here. 


From  (3.1)  and  (3.4) 


*/=  J  P/o- 1  •  •  ■  Pi- 1  |j  -  F.(/.)J/.^/.)  — 

/*7)  J  ~  1 

0  P*P)(l*)e-r°,°F.(U)l.i<tt<,)dU 
JO  A  +  » 

where  f <,('«)  =  1  —  Fa(ta)  and  the  finite  upper  limit  is  set  by  the  fact  that  initiation  l,>  0. 
This  expands  into 

-j: 


P(T]0+i  -to)...  P(T,1  -  r#)F(».)e-^<>[  I  -  Fa(U)]Mto)dU 


(4.11) 


where  the  factors  P,  P  depend  implicitly  on  ta  also  and 

to  =  T)  —  ta. 

Now  we  know  that 

j 

MO  =  2  e~roTk  Rkgt(t  ~  Tk)' 

k  =  0 


go  =  MO, 


(4.11  A) 


whence 


CT, 

Ho 


P(Tjo+i  -  7)  +  /,).. .  P(T/-i  —Tj  +  U)P(t.)Pa{t.) 


Jo 

x  ^  e  roTkRtgt(T,  -  7*  +  ta)  dt.  (4. 1 2) 

This  is  essentially  an  infinite  set  of  recursive  linear  equations  for  Rj  (R0  —  1).  If  their  matrix 
form  is 

R  =  BR  B  =  [bti] 

then  bij  =  I  P(Tj0+i  —  7)  +  to) .  ■  ■  P(Tf_i  —  7)  +  to)P{ta)F i(/«) 


j: 


=  o 


X  e~r°T‘g<(Tj  —  7V  +  ta)dta  if  i  <  j;  f«,  7)  -*j0, 

if  i  >  j. 


(4.I2B) 


all  obtainable  in  principle  by  simple  quadratures.  Since  R0  is  unity  these  equations  are  not 
homogeneous. 


5.  EMBEDDED  MARKOV  PROCESS 

In  another  interpretation  the  by  are  transition  probabilities  (rejection  and  repair)  in  a 
transient  discrete  Markov  process  embedded  in  the  attrition  and  fatigue  process.  The  dement 
bn  represents  probabilities  of  structures  repaired  after  time  7)  being  rejected  again  at  the 
inspection  7). 

However,  in  the  equation  R  =  BR,  the  matrix  B  consists  of  transition  probabilities  among 
an  infinite  number  of  states.  Nevertheless,  these  states  are  still  incomplete,  not  including  attrition, 
hijacking  or  even  new  cracking.  Furthermore,  the  infinite  vector  R  is  described  by  a  single 
transition  B  covering  all  time.  It  is  more  convenient  to  consider  transitions  just  after  each 
inspection  time  Tr,  the  states  at  such  times  may  be  regarded  as  agglomerations  of  the  states 
represented  by  by.  Thus,  structures  cracked  at  7)  include  rejections  and  repairs  from  previous 
inspections. 
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In  this  way,  the  process  R  may  be  identified  with  a  series  of  transitions  of  a  smaller  process 
after  each  inspection.  From  the  Venn  diagram  of  Figure  I ,  the  number  of  states  in  this  smaller 
process  is  six-three  for  cracked  and  three  for  uncracked  structures.  One  of  these  states  has  zero 
probability  (although  it  may  be  considered  an  absorbing  state)  so  that  a  5  x  5  transition  matrix 
is  required  for  each  inspection  time.  These  will  be  denoted  B  for  time  T,.  The  state  ignored  is 
the  rejection  of  uncracked  structures  at  inspection. 

To  find  B*.  it  is  necessary  to  consider  the  possible  combinations  of  various  states.  This 
amounts  to  gathering  previous  results  and  formalising  assumptions,  implicit  or  otherwise.  We 
consider  the  epochs  Ti  +  just  after  any  inspection  and  imagine  the  “mainstream”  of  structures 
as  having  survived  the  initial  hijacking  at  rate  r0,  becoming  cracked  by  time  7i,  and  then  being 
subdivided  into  the  mutually  exclusive  states  of  “immediate"  (i.e.  at  Ti  +  )  rejection,  attrition 
in  the  i-th  interval  (7<,  Ti ,  i]  and  the  mainstream  for  7j,|.  The  last  will  be  augmented,  as  at  each 
inspection,  by  repaired  structures  cracking  again;  one  of  the  state  transitions  to  be  considered. 


5.1  State  Probabilities 

We  begin  by  considering  the  state  probabilities  at  7j  »  which  may  also  be  normalises  for 
the  conditional  transition  probabilities.  Comparison  with  7j«i  +  and  the  use  of  known  results 
then  provide  elements  of  B.  Let  <MO  be  the  integrated  attrition  and  abbreviate  <M.Tt),  F4Tt )  to  <(> 
and  F.i.  In  addition  C,  A  and  R  are  sets  of  respectively  cracked,  failed  or  hijacked  and  rejected 
structures;  before  time  Ti  +  in  the  first  instance.  Set  complements  are  subscripted  c  and  unions 
and  intersections  are  most  conveniently  denoted  by  +  and  the  product  convention. 

Among  the  eight  factors  of  (C  f  Cr)(A  f  AcM.R  f  Rc)  three  are  empty,  namely: 

CAR  Inspection  and  rejection  of  a  cracked  structure  after  attrition; 

CcAR  Inspection  and  rejection  for  hijacked  structures;  and 

CeA'R  Rejection  of  ordinary  uncracked  structures. 

After  obvious  condensations  this  leaves  the  universe 

E  —  CR  f  CA  f  CAcRc  f  CcA  -f  CcAc  (5.1) 

in  which  C,  A  and  their  complements  describe  the  end  result  of  past  history  but  R,  Rc  refer 
to  the  one  inspection,  beginning  the  interval  (Ti,  Tin). 

At  time  7V+  we  know  that 

Pr(CR)  =  Ri  from  (3.3)  with  f  =  Ti  (5.2) 

The  basic  hijack  allowance  states  that  the  mainstream  CrAr  of  uncracked  structures  is 
exp(—  r<,7j)(l  -  F«). 

Without  hijacking,  Pr(CeAr)  would  reduce  to  I  —  F.i  so  that  the  difference 

Pr(CcA)  =  (1  -  f^xi  -  F„),  (5.3) 

the  fraction  for  hijacked  but  uncracked  structures. 

Hence  for  those  already  cracked 

Pr(CR  +  CA  +  CA'Rc)  -=  F.u  (5.4) 

the  cracked  mainstream. 

By  definition  <f>i  =  Pr(A)  =  Pr(CeA  +  CA) 
from  which  subtraction  provides 

Pr(CA)  =  $t—(l  —  e~r°Tl)[l  -  F«]  (5.5) 

Now 

Pr  (cracked  mainstream)  =  Pr(CA  +  CAeRc  +  CR) 
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and  by  subtraction 

Pr{CA'Rt )  -  F.{  R,  t  (I  e  r"T' )fl  F.,)  <P,  1  -e  '«•*■.  (|  F„)  R,  -  0,  (5.6) 

This  is  not  the  mainstream  at  Ti  t  since  it  will  be  affected  by  new  cracking  and  attrition  during 
the  interval  (Tt,  Tt, i). 


FIG  5.  STATES  AFTER  EACH  INSPECTION 


5.2  Transition  Probabilities 

Equations  (5.1)  to  (5.6)  describe  the  five  states  shown  in  Figure  5  and  their  probabilities. 
Now  consider  transitions  between  them  during  (Ti,  7V«i].  The  main  difference  from  the  pre¬ 
viously  used  infinitesimal  Markov  process  is  that  the  finite  inspection  interval  allows  transition 
from  repaired  or  uncrncked  structures  to  any  state. 

Obviously  attrition  is  always  an  absorbing  state.  The  process  as  a  whole  is  therefore  transient 
and  the  ultimate  fate  of  all  structures  is  either  retirement,  or  attrition  by  failure  or  hijacking. 
We  shall  now  examine  transitions  out  of  the  five  basic  conditions  using  set  theory  and  pro¬ 
bability  where  necessary. 

The  more  difficult  elements  will  require  interpretation  of  various  terms  in  the  expressions 
(3.5),  (3.6)  or  (3.3)  for  attrition  or  rejection.  Integrals  of  these  may  appear  as  convolution 
functions. 


5.11  CR 

in  terms  of  time  variates  this  has  the  total  measure 


m 


PkP(ta)?  r"l“[  I  —  F a{t  Q)\j\(t  0)dt  a ,  to  —  Tt  —  /a. 


which  leads  to  the  component  of  initiation  density  Rtf  r°tto  Tl)gtC"  r'). 

The  arguments  of  Section  3  apply  equally  to  all  components  so  that  in  the  interval  (T(  r<  +  j) 
this  corresponds  to  the  partial  density  of  crack  length 

(RtlR(a))e  r«‘l»  T'’(l  -  F.IMkrt'  T,  /.),  t,  =  t  -  r„  (5.7) 

which  forms  part  of  (3.5). 

Corresponding  to  (3.6)  one  finds  the  attrition  rate 

^(/|From  /(<)  =  Rtf r»,f  r'>[l  GiO  —  7V)J 


f/  -  T, 

f  Rt  I  e  '»“«  Ti'gtd  T,-t.)dFJl.) 
JO 
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from  which  integration  shows 


RiPr(CRt  ->  CeA\Ri)  = 


fr<+ 1 


rae  toit-Tn[\  _  c,(/  -  T,)\dt 


(5.8) 


and 


RiPr(CRt  ->  Ot!*()  = 


nt  -  T, 
0 


e  ^Uo-T,)gi{l  _  r,  -  r.)  dFJ,u)dt  (5.9) 


The  first  of  these  may  be  integrated  by  repeating  the  general  argument  for  (5.3).  The  transition 
element  is  therefore 


Pr(CRt  -*  CcAlRi)  =  (I  -  e  '^Ti)  (l  _  G,(ATt)),  (5.10) 

where  AT,  =  7m  —  Tt. 

Equation  (5.9)  has  no  simpler  form  but  it  may  be  interpreted  as  a  probability.  Since 

1  -  Hi(t)  =  e~r°*[\  -  G,(r)]  (S.10A) 

is  a  survival  probability  under  attrition,  its  complement  is  a  true  distribution  function  and  never 
defective.  Then,  as  a  convolution, 

Pr(CR,  -+  .41*0  =  Pr(t0  +  /«  <  AT,) 

=  H,*F.(AT,),  say, 

where  /«  ~  Ht(t0),  Hi  from  above,  and  ~  Fa(ta) 

for  the  reliability  based  crack  growth  time.  This  includes  R,  -*>  C,A  but  the  time  does  not  allow 
inspection.  By  subtraction  of  (5.10) 

Pr(CR,^CA\R,)  =  H,*F.(AT,)-0  -e  ^Tl)(l  _  G,(AT,)),  (5.11) 

we  have  already  noted  that 

Pr(CR,  -*  CcAe\Ri)  =  e~roJT‘(l  -  G,(AT,))  (5.12) 

augmenting  the  mainstream  of  uncracked  structures. 

There  remain  two  more  transitions  to  CR,+\  or  CAcRe,  in  the  cracked  mainstream.  The 
first  of  these  follows  in  the  same  way  as  Rt  itself,  but  from  the  partial  density  (5.7).  Thus 

CAT, 

Pr(CRi  ->  CR\R,)  =  e  'o^-TOp(ta)[l  -  Fa(t.)]g,(AT,  -  t.)du 

Jo 

If  we  recall  our  assumption  P(0)  =  0,  this  takes  the  form 

Cat, 

P(ta) P  -  F.(u)]dH,(Ar,  -  ta)  =  (PU  -  F.])*H,(AT,)  (5.13) 

Jo 

Finally,  for  transitions  to  the  cracked  mainstream, 

Pr(CR,  -*  CAcRAR,)  =  1  -  Eq.  (5.10)  -  ...  -  (5.13) 

=  H,(AT,)  -  -  P)F.  +  PXAT,)  (5.14) 

=  ((1  -  P)(l  -  Fal)*Ht(AT,). 


5.2.2  CAcRc  =  M 

This  is  defined  by  (5.6)  for  entry  to  the  inspection  interval  (7«,  7ml.  Just  after  7m,  it  has 
changed  by  rejection  and  an  influx  of  structures  repaired  after  the  previous  inspection  T,. 

Obviously  M  -*■  CeA  and  M  -*  CcAc  are  impossible;  also  the  influx  of  cracked  structures, 
newly  repaired  and  otherwise,  complicates  Pr(M,  -*■  M,+\\M,). 
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Equation  (S.6)  provides  the  unconditional  Pr(Mt+i)  as  well  as  Pr(Mt)  but  this  influx  makes 
it  more  convenient  to  start  again  from  the  crack  length  density  (3.4).  At  time  7<+j,  the  only 
cracks  from  Mi  must  be  those  for  which  U  >  AT{.  If  (3.4)  is  restricted  to  these  at  t  =  7Wi, 
then  all  attrition  during  (7V,  7<  +  i]  is  accounted  for,  and  also  the  rejections  Ri+i.  Then,  changing 
the  cracked  length  variate  to  /«, 


Pr(Mt)Pr(M,+i\M,) 


(I  -  F*U))e-'*> 


~T,)dta 


with  to  =  Tt+i  —  to.  This  differs  from  Pr(Mi+i)  only  by 


cat, 

N,=  I  (1  -  /*(/.))[!  -  Faito)\e  ^of*(Tul  -  ,.)dt.  (5.15) 

Jo 

after  some  obvious  simplifications.  As  in  (5. 1 1)  it  would  now  be  possible  to  define  a  distribution 
1  —  (exp  —  r0/Xl  —  F4.I)),  and  thus  define  this  term  as  another  convolution.  However,  this 
includes  elements  of  CeA,  the  attrition  of  uncracked  structures,  and  to  abbreviate  Ni  we  use 
the  defective  density  hJ^to)  —  (exp  —  r0t0)f4to)  so  that 

Ni  =  ((I  —  P)[\  —  Fo])»htAT,)  (5. 1 5A) 


If  we  are  guided  by  context  we  may  conveniently  let  Pr(Mt)  =  Mi,  using  the  same  symbols 
for  the  event  and  its  probability.  Then,  using  (5.6)  and  (5.15), 


Pr(Mt  -+  Mi+\\Mi)  =  ( Mi+l  -  Nt)IMi  = 


1  —  e  f»T'+'ll  —  F*+i]  —  Ri+t  —  —  Ni 


(516) 


The  difference  Mt  +  Nt  —  Mi+i  must  now  be  divided  among  Mt  -*■  CRm  and  Mi  -*  CAu  i 
to  find  the  corresponding  transition  probabilities.  Consider  the  attrition  of  cracked  mainstream 
structures  during  (7V,  Ti+i]  given  by  (3.6),  and  based  on  (3.5).  In  the  latter,  CAcRc  or  M( 
correspond  to  use  of  the  partial  density  f«\(ta)  during  the  i-th  inspection  interval.  This  excludes 
attrition  of  structures  repaired  at  Ti.  If  they  are  cracked  however,  their  absolute  probability 
from  (5.1 1)  is  known  to  be 

RiH(*F.(AT,)  -  (5.10)  =  RtPr(CRt  -+  CA\R,). 

In  (3.6),  the  second  term  is  the  attrition  density  of  previously  cracked  structures,  and  those 
repaired  at  7).  Writing  the  second  term  as 

<£.  -  fo(exp  -  /w)[l  -  F^t)\ 

and  allowing  for  (5.16) 

rr(+ 1 

M,Pr(CAeRc  ->  CA\Tt  <  /  <  Tt+i)  =  I  (<M0  -  rtf ^[1  -  FJ(t)])dt  + 

+  Ri( I  -  e  'oJT/fi i  _  Ci(ATi)]  -  R,Ht*F.(AT,) 

The  second  term  of  the  integral  follows  as  before  from  the  increment  of  (5.4).  After  some 
reduction  this  becomes 

M(Pr(M  -*  C4|/-th)  =  A0t  +  AF«( I  -  e-*Tf* ')  -  RtH,*F^AT,) 

-  (I  -  e  ,oAT>){e-T»Tt[\  -  F„]  -  Ri[\  -  Gt(ATt)])  (5.17) 

From  the  remarks  below  (5.15) 

MtPr(Mi  -*  CffmlMh)  =  M(  +  Ni  -  M,+l  -  Eq.(5.17) 
which  eventually  reduces  to  MtPr(Mi  -*  C/?i+i|/-th)  = 

-Rt~  AF«  +  Ni  +  R,{H,*F^AT,)  -  (I  -  e^OCl  —  C7«(A7*)P 
with  Ri,  Ni  and  Hi  defined  by  (3.1)  or  (4.1 1)  -  (4.12),  (5.15)  and  (5.10A). 
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5.2.3  C,A, 


This  is  the  only  remaining  non-trivial  state.  Any  transition  is  possible  from  it  and  its  absolute 
probability  at  Tt  I  is 

Pr(CeAc)  e  r°T,(l  -  F.<] 

During  (7),  ,  1  ]  Pr(CrA)  increases  by  the  increment  of  (5.3)  and  this  comes  from 

CeAr  t  CRi.  Thus, 

e  '“r'[l  -  F.f]Pr(CcAe  -  GAii-th)  -  APr(CcA)  -  R,Pr(CR,  -  CcA\R,) 

—  (!-«*  r»,T'){e  '“T'll  -  F„]  -  K([l  -  CddT,)]}  -  AF«( I  -  i)  (5.19) 
using  (5.3)  and  (5.10). 

As  a  whole,  G  comes  from  CcA(  |  CRi  and  of  course  its  increment  has  measure  AF.t,  the 
decrease  in  G-  For  the  increment  from  CeAc,  subtract  the  transitions  from  C7t<  giving 

Pr(CtAc  -  G»i)  -  AF.i  -  R,Pr(R,  .  C(A  I  R  \  M)\Ri) 

-  AF.t  -  Ri{H,{ATi)  -  (I  e  ^T<)t»  -  G(dr()]}  (5 .20) 

which  follows  from  (5.1 1),  (5.13)  and  (5.14)  respectively. 

Consider 

-  APr{CcAc)  =  Pr(CcAc)\Pr{CeA<  -  C,u 7-th)  -f  Pr(CcAc  -*  CA.A  +  /(,);/- th) 
whence  by  transposing  (5.19)  and  (5.20) 

c  'oT,[ I  _  F.t)Pr(C'Ac  '  CcAc  i-th)  =  RtH,(AT,)  AF.,e  '°'T‘(  1  -  e  '°T>)  (5.21) 

with  Hi  from  (5.I0A). 

As  was  done  with  transitions  CRi  *  G  i  we  must  now  partition  (5.20).  After  subtracting 
the  correction  term  in  (5.19)  for  Rt  *  G.  i,  one  may  imagine  the  three  types  of  transition 
during  (7),  7*.i]  being  driven  by  /.<  i.  Attrition  is  simplest  with  a  contribution  from  part  of 
(3.6)  minus  (5.1 1).  Thus,  the  appropriate  rate  is 

ft-  Tt 

<k<  Ti  /  r,.,)  =  e  '»>'  -  ta)dFaUa ) 

Jo 

I  T, 

e  -  RtgiV o  -  Ti))dFa(u) 

0 

Over  the  inspection  interval  this  integrates  to 
e  roT,[ i  _  F.i]Pr(CcA'  ►  CA  i- th) 

=  A<Pi  -  e  r“Tl{ I  _  r„ir,)(|  /r.,j  _  Rt{Ht* Fa(ATt)  -  (I  -  er°'T‘) [I  -  Gi(ATt)]}  (5.22) 

using  (5.1 1)  for  the  subtracted  term. 

Inspection  of  (3.5)  for  /  -  7"(,i  f  indicates  the  crack  length  density 

R(a]f(a\T„i)  =  e  '<*>[1  -  F.(r„)KI  -  P{t.))f.i  i(to) 

for  structures  from  CrAc.  For  rejected  structures,  we  replace  I  —  P  by  its  complement  and 
integrate  to  obtain 

e (I  -  F.i]Pr(C'Ac  >  CR,>i\i-lh)  -  </*(!  -  F.))*H^ATt)  -  *,(/»[ I  -  F.])*H,(ATt)  (5.23) 

using  (5.13)  with  the  assumption  P(0)  =  0  and  (5.I5A),  (5.I0A)  for  definitions  of  H.  and  Ht. 

We  now  know  probabilities  for  transition  to  CA  f  Crtr+i  and  to  G+t  as  a  whole  from 
equations  (5.22),  (5.23)  and  (5.20).  Since  the  transition  events  are  disjoint 

e  ,oTi(l  —  FH)Pr(C'Ac  *  A/m;/-th)  =  AF.i  —  Aty  ■+  e  r®T<(I  —  e  r»3T,)[l  —  fil 

-  RiNi  -  (/II  -  F.])*H4AT,)  (5.24) 

using  (5.  ISA). 
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5.3  Tnumitioa  Matrices 


It  is  now  possible  to  arrange  the  elements  above  into  transition  matrices  for  each  epoch 
Ti  + .  In  many  cases,  normalising  to  conditional  probabilities  introduce  awkward  fractions  so 
that  it  is  most  convenient  to  place  the  normalises  in  a  diagonal  prefactor.  Doing  this,  and 
gathering  elements  from  above,  leads  to  the  matrices  shown  in  Tables  2  and  3. 

Let  P  1  =  [/*<]■ 1  be  the  prefactor  where  only  Pa  and  Pi  differ  from  unity. 

From  Section  S.l  we  know  the  absolute  state  probabilities  of  Table  2  which  provide 
Pa  and  Pi. 


TABLE  2 
Normalising  Factors 


j 

State 

Name 

Absolute  Probability 

D 

1 

CcA 

Hijacked 

£ 

1 

i-' 

1 

i 

2 

CA 

Failed  or  hijacked 

0(  -  (1  -  e*r"r()[l  -  F.<] 

i 

3 

CR 

Rejected 

Ri 

i 

4 

CA'R'  =  M 

Cracked  mainstream 

1  e  r°Ti[\  -  F„) 

n 

5 

CcAc 

Uncracked  mainstream 

e~r°Ti[l  -  F.«] 

H 

*  As  in  previous  column. 


In  previous  derivations,  all  results  depend  on  functions  which  are  fully  defined  at  time 
Tt+  on  the  range  of  crack  lengths  Ou{fl0,  oo)  which  imply  the  other  functions  needed  from 
(T,  Ti+i].  This  validates  the  Markovian  nature  of  the  discrete  transitions. 


6.  CONCLUSIONS 

Using  the  same  procedures,  the  previous  analysis3  for  fatigue  life  distributions  in  one-crack 
models  with  hijacking  have  been  extended  to  cases  with  arbitrary  renewals  at  general  inspection 
intervals. 

The  conclusions  are  similar  to  those  derived  previously  if  one  postulates  a  defective  distri¬ 
bution  /.  of  crack  initiation  lives  as  affected  by  renewals,  and  includes  factors  depending  on 
the  operating  characteristic  of  the  inspections.  The  moment  generating  function  has  the  same 
form  as  before3  and  with  f0(t0)  replaced  by  /.(/<>)•  The  proof  of  this  is  facilitated  by  treating 
inspection  and  crack  life  density  together. 

The  two-stage  fatigue  process  is  Markovian  in  continuous  time.  At  any  time  a  structure 
must  be  in  one  of  five  states  corresponding  to  combinations  of  attrition,  cracking  and  rejection 
at  inspection  or  their  complements.  Transition  between  these  states  is  a  variable  Markov  chain 
embedded  in  the  continuous  process  with  epochs  conveniently  placed  just  after  inspections. 

6.1  Implementation 

Life  distributions,  risk  rates  and  rejection  probabilities  may  be  found  from  two  FORTRAN 
IV  programs  developed  from  the  preceding  theory.  The  second  of  these  has  options  for  random 
crack  rates  and/or  run-time  setting  of  inspections.  This  generalisation  and  the  program  will  be 
described  in  two  further  reports. 
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Harvard  Library,  Engineering  217 

Illinois  Professor  N.  M.  Newmark  218 

Iowa  State  Dr  G.  K.  Sekory,  Mechanical  Engineering  219 

Johns  Hopkins  Department  of  Mechanical  Engineering  220 

Massachusetts  Library,  School  of  Engineering  221 

Institute  of 

Technology  Library  222 

New  Mexico  Dept  of  Nuclear  Engineering,  Professor  J.  A.  Horak  223 

Pennsylvania  Library,  School  of  Engineering  224 

Stanford  Library,  Department  of  Aeronautics  225 

Wisconsin  Memorial  Library,  Serials  Department  226 

Spares  227-236 


